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(54) A method of using inductance for determining the position of an armature in an 
electromagnetic solenoid 

(57) An improved method for controlling the landing 
velocity of an armature in an electromechanical actua- 
tor, such as a fuel injector, fuel pressure regulator, or 
engine valve actuator is provided. The position and 
velocity of an armature during a stroke is dynamically 
estimated by calculating the inductance and rate of 
change of inductance of the actuator coil in real-time as 
the armature moves through its stroke, compensating 
for non-linear permeability and magnetization effects 
due to changing gap, temperature, magnetic material 
properties or magnetic architecture, normalizing the cal- 
culated inductance value at the end of a stroke (zero 
gap), and mapping the value of normalized inductance 
to correspond to an armature position by an algebraic 
transformation. Inductance may be used directly as a 
position variable without mapping it to units of position. 
Rate of change of inductance may be used as a rate 
variable without mapping it to units of velocity. 
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Description 

[0001] This application claims the benefit of U.S. 
Provisional Application No. 60/143,619, filed July 13, 
1999, which is hereby incorporated by reference in its 

entirety. 

Background of the Invention 

[0002] This invention relates to a high-speed, high- 
force electromagnetic actuator, and particularly to an 
electromagnetic actuator and method for opening and 
closing a valve of an internal combustion engine, driving 
a high pressure fuel injector, or operating a high pres- 
sure fuel regulator. More particularly, this disclosure 
relates to an apparatus and method of dynamically 
measuring the inductance and rate of change of induct- 
ance of a, electromechanical actuator as the armature 
moves from one pole piece toward another and inferring 
armature position and velocity from the measured 
inductance. Still more particularly, this invention relates 
to an electronic apparatus and method of using induct- 
ance and rate of change of inductance for dynamically 
controlling the landing velocity of an armature in a fuel 
injector or an electromagnetic actuator for opening and 
closing a valve of an internal combustion engine. 
[0003] Electromagnetic actuators, such as fuel 
injectors, actuators for opening and closing a valve in an 
internal combustion engine (hereinafter "Electronic 
Valve Timing" or "EVT* actuators), and fuel pressure 
regulators, typically include a solenoid for generating 
magnetic force. A solenoid is an insulated conducting 
wire wound to fomi a tight helical coil. When current 
passes through the wire, a magnetic field is generated 
within the coil in a direction parallel to the axis of the 
coil. The resulting magnetic field exerts a force on a 
moveable ferromagnetic armature located within the 
coil, thereby causing the armature to move from a first 
position to a second position in opposition to a force 
generated by a return spring. The force exerted on the 
armature is proportional to the strength of the magnetic 
field and the strength of the magnetic field depends on 
the number of turns of the coil and the amount of current 
passing through the coil. 

[0004] While it will be appreciated by those sl<i!led 
in the art of electromechanical actuators that the tech- 
niques described in the present disclosure may be 
applied to any electromechanical actuator, including, for 
example, fuel injectors or fuel pressure regulators, for 
purposes of clarity the present invention will be 
described primarily in the context of an EVT actuator for 
opening and closing a valve of an internal combustion 
engine. 

[0005] An EVT actuator generally includes an elec- 
tromagnet for producing an electromagnetic force on an 
armature. The armature is typically neutrally-biased by 
opposing first and second return springs and coaxially 
coupled with a cylinder valve stem of an engine. In oper- 



ation, the armature is held by the electromagnet In a first 
operating position against a stator core of the actuator. 
By selectively de-energizing the electromagnet, the 
armature may begin movement towards a second oper- 

5 ating position under the influence of a force exerted by 
the first return spring. Power to a coll of the actuator 
may then be applied to move the armature across a gap 
and begin compressing the second return spring. 
[0006] As can be appreciated by those skilled in the 

10 art, it is desirable to closely balance the spring force on 
the amnature with the magnetic forces acting on the 
armature in the region near the stator core so as to 
achieve a near-zero velocity "soft landing" of the arma- 
ture against the stator core. In order to obtain a soft- 

15 landing of the armature against the stator core, power to 
the coil may be modulated to reduce the armature 
velocity as the armature approaches the stator core in 
the second position. The coil may then be re-energized, 
just before landing the armature, to draw and hold the 

20 armature against the stator core. In practice, a soft land- 
ing may be difficult to achieve because the system is 
continually perturbed by transient variations in friction, 
supply voltage, exhaust back pressure, armature center 
point valve lash, engine vibration, oil viscosity, toler- 

25 ance stack up, temperature, etc. 

[0007] Soft landing techniques are becoming espe- 
cially important with modern high-pressure fuel injectors 
and direct injection fuel injectors that employ strong 
return springs. Soft landing the injector amriature 

30 reduces injector noise and internal wear. In addition to 
noise reduction, soft landing has the benefit of reducing 
power consumption in the actuator because it enables 
controlled metering of the coil current so as to only 
place the required amount of magnetic energy in the 

35 system necessary to actuate the armature. Soft landing 
techniques may also be applied to control the landing 
velocity of an armature in a high pressure fuel regulator. 
[0008] In the case of EVT actuators, experimental 
results for particular engines and actuator arrange- 

40 ments indicate that to achieve quiet EVT actuator oper- 
ation and prevent excessive impact wear on the 
armature and stator core, the landing velocity of the 
armature should be less than 0.04 meters per second at 
600 engine rpm and less than 0.4 meters per second 

45 6,000 engine rpm. In order to achieve these results 
under non-Ideal conditions (e.g., the harsh environment 
of an internal-combustion engine), it is necessary to 
dynamically monitor and adjust the magnetic flux gener- 
ated within the magnetic circuit to compensate for van- 
so ations in operating voltage, friction within the actuator, 
engine back-pressure and vibration, during every stroke 
of the amriature. External sensors, such as Hall sensors, 
have been used to measure flux in electromagnetic 
actuators. However, sensors have proven to be too 

55 costly and cumbersome for practical applications. 

[0009] FID (proportional, integral, derivative) con- 
trol methods have been proposed to control the landing 
velocity of an armature in an electromagnetic actuator. 
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An example of using PID methods to control the landing 
velocity of an armature in an electromagnetic actuator is 
disclosed in U.S. Patent Application No. 09/434.513. 
filed November 5. 1999 and entitled "Method of Conrv 
pensation for Flux Control of an Electromechanical 
Actuator." the contents of which is hereby incorporated 
in its entirety into the present specification by reference. 
Generally, PID control systenns can only perfectly com- 
pensate a linear system with state variables that are not 
interactive. Electromagnetic actuator systems are, how- 
ever, highly non-linear due at least in part to changing 
magnetic permeability as the armature moves within the 
actuator In addition, the state variables of an actuator 
(i.e., flux, position, and velocity) are highly interactive. In 
order to apply PID methods to control the landing veloc- 
ity of an armature in an electromagnetic actuator, sim- 
plifying linear approximations are necessary, e.g., the 
system must be presumed linear over small armature 
displacements and the state variables must be pre- 
sumed to be independent. Accordingly there is a need 
for a true multivariate control system capable of control- 
ling all state variables simultaneously and compensat- 
ing a nonlinear feedback control system. 
[0010] The present invention overcomes the two 
classical limitations of pure PID control described above 
by providing a sensorless position estimator that ena- 
bles automatic calibration of the system. Sensorless 
position estimation accounts for much of the non-linear- 
ity of the system. Knowing armature position throughout 
the armature stroke makes it possible to self-calibrate 
the control system. This is because once amnature posi- 
tion is known, together with another state variable such 
as velocity, it is possible to employ known non-linear 
multivariate feedback control algorithms to control the 
system. 

[001 1 ] The prior art lacks a practical and cost effec- 
tive method of dynamically measuring armature posi- 
tion during the armature stroke. While lasers have been 
used in laboratory settings to measure armature posi- 
tion, it is not practical or cost effective to put a laser on 
actuators manufactured for large-scale production. 
Other more cost-effective methods of position sensing 
have not proven to be accurate and durable enough. For 
example, in automotive applications, position sensors 
must be able to withstand the temperature and vibration 
extremes of being mounted on an engine. Sensor- 
based techniques also present the problem of cabling 
the signal through a potentially electrically noisy envi- 
ronment. Accordingly, there is a need to estimate arma- 
ture position in a sensorless manner. 
[001 2] Thus, a need exists for a sensorless self-cal- 
ibrating control system and method for an electromag- 
netic actuator capable of dynamically compensating for 
non-ideal disturbances that exist in and near internal 
combustion engines. Further, a need exists for a high- 
speed sensorless control system and method for an 
electromagnetic actuator capable of detecting and com- 
pensating for the above-described non-ideal conditions 



during each stroke cycle of the armature. 
Summary of the Invention 

5 [0013] A sensorless method of controlling the land- 
ing velocity of an armature in an electromagnetic actua- 
tor is provided. The method disclosed dynamically 
measures actuator inductance and rate of change of 
inductance as the armature moves within the coil. The 

10 B-H magnetization characteristics of the actuator during 
an armature stroke are determined during actuator 
operation and the measured inductance and rate of 
change of inductance are thereby compensated for non- 
linear permeability and magnetization effects. The 

75 measured inductance may be normalized at zero gap. 
In a preferred embodiment, the normalization at zero 
gap is to unity (1 .0). From inductance, an estimation of 
position is made; from rate of change of inductance, 
armature velocity information is inferred. The armature 

20 position and rate information are provided to a control 
system for modulating a current delivered to the actua- 
tor, thereby controlling the armature landing velocity. 

Brief Description of t/ie Drawings 

25 

[0014] The accompanying drawings, which are 
incorporated herein and constitute part of this specifica- 
tion, illustrate presently preferred embodiments of the' 
invention, and, together with the general description 
30 given above and the detailed description given below, 
serve to explain features of the invention. 

Rgure la illustrates a sectional view of an elec- 
tronic valve timing electromagnetic actuator pro- 

35 vided in accordance with the principles of the 
present invention, shown in a valve open position. 
Rgure lb illustrates a sectional view of an elec- 
tronic valve timing electromagnetic actuator pro- 
vided in accordance with the principles of the 

40 present invention, shown in a valve closed position. 

Rgure 2 illustrates a sectional view of a direct injec- 
tion fuel injector provided in accordance with the 
principles of the present invention. 
Rgure 3 is a system block diagram in accordance 

45 with a preferred embodiment of the present inven- 
tion. 

Rgure 4 illustrates the relationship between coil 
voltage and magnetic flux density in accordance 
with a preferred embodiment of the present inven- 

50 tion. 

Rgure 5 is a schematic diagram illustrating a 
method of dynamically determining the inductance 
of an electromagnetic actuator as the armature 
moves from one pole piece to another, in accord- 

55 ance with a preferred embodiment of the present 
invention. 

Rgure 6 illustrates the waveforms representing 
measured coil current and voltage, and calculated 
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coil inductance using digital signal processing tech- 
niques in accordance with a preferred embodiment 
of the present invention. 

Figure 7 illustrates typical B-H magnetization 
curves over a range of air gaps. 
Figure 8 illustrates mu factor autocalibration in 
accordance with the present invention. Figure 9 
illustrates the results of sensorless armature posi- 
tion estimation in accordance with the present 
Invention. 

Figure 10 illustrates a comparison of inductance 
with the integral of magnetic flux. 
Figure 11 illustrates normalized inductance and 
rate of change of inductance determined in a sen- 
sorless manner in accordance with the present 
invention. 

Figure 12 is a block diagram of a lookup table 
implementation for determining running set points. 
Figure 13 illustrates a comparison of ideal induct- 
ance with measured inductance in accordance with 
the present invention. 

Figure 14 illustrates a sensorless soft landing of an 
amnature in an electromagnetic actuator in accord- 
ance with the present invention. 

Detailed Description of ttie Preferred Embodi- 
ment(s) 

[0015] The present invention will be described pri- 
marily in relation to an EVT actuator. However, as will be 
appreciated by those skilled in the art, the present 
invention is not so limited and may be applied to any 
type of electromechanical actuator including, for exam- 
ple, fuel injectors and fuel pressure regulators. 
[0016] In accordance with a preferred EVT embodi- 
ment, Figures la and lb illustrate an electromagnetic 
actuator 10 for opening and closing a valve in an inter- 
nal combustion engine. The electromagnetic actuator 
10 includes a first electromagnet 12 that includes a sta- 
tor core 14 and a solenoid coil 16 associated with the 
statorcore 14. A second electromagnet 18 is disposed 
in opposing relation to the first electromagnet 12. The 
second electromagnet includes a stator core 20 and a 
solenoid coil 22 associated with the stator core 20. The 
electromagnetic actuator 10 includes an armature 24 
that is attached to a stem 26 of a cylinder valve 28 
through a hydraulic valve adjuster 27. The armature 24 
is disposed between the electromagnets 12 and 18 so 
as to be acted upon by the electromagnetic force cre- 
ated by the electromagnets. In a de -energized state of 
the electromagnets 12 and 18, the armature 24 is main- 
tained in a neutrally-biased rest position between the 
two electromagnets 12 and 18 by opposing return 
springs 30 and 32. In a valve closed position (Figure 
lb), the armature 24 engages the stator core 14 of the 
first electromagnet 1 2. 

[0017] To initiate motion of the armature 24 and 
thus the valve 28 from the closed position into an open 



position (Figures 1a & 1b), a holding current through 
solenoid coil 16 of the first electromagnet 12 is 
removed. As a result, a holding force of the electromag- 
net 12 falls below the spring force of the return spring 30 

5 and thus the armature 24 begins moving under the force 
exerted by return spring 30. It is necessary to build 
enough magnetic flux in the coil 22 so there will be suf- 
ficient magnetic force to make the armature 24 move 
from one stator 14 to another 18 while overcoming the 

10 opposing neutrally-biased return springs. To catch the 
armature 24 in the open position, a catch current is 
applied to the electromagnet 18. Once the armature has 
landed at the stator core 20, the catch current is 
changed to a hold current which is sufficient to hold the 

75 armature at the stator core 20 for a predetermined 
period of time. It is desirable to dynamically control the 
catch current to achieve a near-zero velocity "soft" land- 
ing of the armature against the stator core. 
[001 8] An example of using rate of change of flux as 

20 a feedback variable is taught in U. S. Patent Application 
No. 09/025,986. filed February 19, 1998 and entitled 
"Electronically Controlling the Landing of an Armature in 
an Electromagnetic Actuator", the contents of which is 
hereby incorporated in its entirety into the present spec- 
ks ification by reference. 

[0019] An example of feedback control based on a 
rate of change of flux without the need for a flux sensor 
is disclosed in U. S. Patent Application No. 09/122,042, 
filed July 24, 1 998 and entitled "A Method for Controlling 

30 Velocity of an Armature of an Electromagnetic Actua- 
tor," the contents of which is hereby incorporated in its 
entirety into the present specification by reference. 
[0020] According to a presently preferred embodi- 
ment, an improved apparatus and method for controlling 

35 the landing velocity of an amnature in an electrome- 
chanical solenoid, such as an EVT actuator or a fuel 
injector will now be described. Referring to Figures 1-3, 
the position of the armature 24 during a stroke may be 
dynamically estimated by calculating the inductance of 

40 the actuator solenoid in real-time as the armature 24 
moves through its stroke; compensating for non-linear 
permeability and magnetization effects due to changing 
gap; normalizing the calculated inductance value to 
always equal unity (1.0) at the end of a stroke (zero 

45 gap); and mapping the value of normalized inductance 
to correspond to an armature position by an algebraic 
transformation. In a preferred embodiment, the induct- 
ance may be used directly as a position variable without 
mapping it to units of position, thus simplifying the 

50 implementation. 

[0021] In similar fashion, the velocity of the arma- 
ture 24 during a stroke may be dynamically estimated 
by calculating the rate of change of inductance of the 
actuator solenoid in real-time as the armature 24 moves 

55 through its stroke; compensating for non-linear permea- 
bility and magnetization effects due to changing gap; 
and mapping the value of rate of change of inductance 
to correspond to armature velocity by an algebraic 
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transformation. In a preferred embodiment, the rate of 
change of inductance may be used directly as a rate 
variable without mapping it to units of velocity, thus sim- 
plifying the implementation. 

[0022] The control loop logic that modulates the coil 
current, and ultimately controls the annature velocity, 
requires as inputs armature position, armature velocity 
and magnetic flux density. Accordingly, in a preferred 
embodiment, armature position may estimated as being 
proportional to a normalized value of inductance and 
armature velocity may be estimated as being propor- 
tional to the rate of change of inductance. 

Dynamic Calculation of Inductance 

[0023] Referring to Figure 4, application of Kir- 
choff's voltage law around the loop yields the following 
relationship: 

V(t) = NdO/dt + l(t)R coil. Equationi 

(where N is the number of turns of the coil, d<l>/dt is 
the rate of change of magnetic flux. I is coil current, 
and Rcoii is not constant). 

[0024] According to a presently preferred embodi- 
ment, a complete processing of the above equation is 
dynamically performed in iterative fashion during actua- 
tor operation. The simplifying approximations of linear- 
ity, independence of state variables (position, velocity, 
and flux density) and the negligible effect of the I R term, 
that were necessary to enable the prior art PID-type 
control, are not necessary in a presently prefen-ed 
approach. In a presently preferred approach, all terms 
of Equation 1 are included in each iterative calculation. 
[0025] In a presently preferred embodiment, com- 
pensation may be made for changes in coil resistance 
due to temperature variations. For example, real time 
resistance measurements may be obtained at the end 
of each armature stroke cycle by measuring the coil 
voltage necessary to maintain a steady-state current 
through the coil and applying Ohm's law to calculate 
resistance. This method of dynamically measuring coil 
resistance is particularly convenient because when a 
steady-state current is applied at the end of an armature 
stroke, d<t>/dt is zero and the voltage drop across the coil 
is IR. With V and I known, R may be readily computed. 
The updated value of R may then be used during the 
next iterative calculation of Equation 1. 
[0026] The basic relationships between magnetic 
flux, <t> rate of change of magnetic flux, d<l>/dt, and 
inductance. L. are as follows: 

<l> = ld<t>/dt, Equation2 
(where <I> is magnetic flux); and 



L = <t>/l, Equations 

(where L is the inductance of the actuator and I is 
coil current). 

5 

[0027] The resistance of the coil may be dynami- 
cally measured during the operation of the electromag- 
netic actuator as follows. The coil voltage may be 
detemnined either by direct measurement or from the 

JO flux mirror circuit method disclosed In U.S. Patent No. 
5,991,143, entitled "Method for Controlling Velocity of 
an Armature of an Electromagnetic Actuator," which is 
hereby incorporated into the present specification by 
reference in its entirety. When a known steady state cur- 

15 rent is applied, the resistance of the coil may be deter- 
mined by applying Ohm's law: R coin^ coi/' steadystate • 
By this method, the resistance of the coil may be 
dynamically measured during each armature stroke. 
[0028] Referring to Figures 3 and 5, the inductance 

20 of the actuator may be dynamically calculated as the 
armature moves from one pole piece to another by solv- 
ing equations 1-3 above in iterative fashion during actu- 
ator operation. . With reference to Figure 5, the 
inductance of the actuator may be computed as follows. 

25 The coil resistance input 52 and coil voltage 50 are 
inputs to the system and may be detemnined by any 
convenient method, including direct measurement or by 
use of the flux min-or circuit method described above. 
As will be appreciated by those skilled in the art, the 

30 direct measurement method requires apparatus suffi- 
cient to detect a small differential voltage in the pres- 
ence of a large common mode voltage, accordingly the 
flux mirror method is preferred. The coil cun-ent 54 is a 
readily measured input to the system because coil cur- 

35 rent 54 is under servo control via a controlled current 
source (not shown). 

[0029] A microprocessor for computing inductance 
L in a dynamic fashion, as described above, must be 
capable of handling a complete cyde of processing and 

40 output the control signal in approximately 40 microsec- 
onds for an EVT actuator, assuming an armature flight 
time of approximately four milliseconds. Figure 5 is a 
schematic representation of a method of computing the 
inductance of an actuator using a commercially availa- 

45 ble microprocessor. An exemplary suitable microproc- 
essor in accordance with a preferred embodiment is a 
TMS320 C3x/4x Digital Signal Processor chip available 
from Texas Instalments. With currently available tech- 
nology, the entire process could feasibly be imple- 

50 mented with many alternative DSP microprocessors, 
digital integrated circuits, or analog integrated circuits. 
In the case of soft landing a fuel injector annature. the 
flight time may be. for example, on the order of 200 
microseconds. Accondingly, with fuel injectors, a high- 

55 speed analog controller is a preferred embodiment in 
order to achieve the necessary processing speed. In the 
case of fuel pressure regulators, a DSP processor may 
be used in a preferred embodiment 
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[0030] Referring again to Figure 5, and in accord- 
ance with equations 1 -3 above, the resistance input 52 
is multiplied 56 by the current input 54, yielding IR, as 
shown symbolically at 58. The calculated value of IR is 
subtracted 60 from the coil voltage input 50, yielding 
rate of change of magnetic flux, d<t>/dt, as shown sym- 
bolically at 62. The flux, 4> 66, is computed by integrat- 
ing the rate of change of flux d<J>/dt 62, as indicated at 
64. Inductance L 70 of the actuator is computed by 
dividing the flux <I> 66 by the coil current input 54, as 
indicated at 68. The Inductance L 70 is then seated 72 
to units of millihenrys (mH). 

Air Gap and Permeability Compensation 

[0031 ] As the armature moves within the solenoid, 
the inductance changes because the reluctance of the 
magnetic circuit is changing due to the changing perme- 
ability of the magnetic circuit. Reluctance in a magnetic 
circuit is analogous to resistance in an electric circuit. 
The components of reluctance are analogous to series 
resistors, a first being of low resistance and correspond- 
ing to the permeability of the ferromagnetic core (arma- 
ture), and a second being of high resistance and 
corresponding to the permeability of air. As the arma- 
ture moves toward the stator core, the total air gap con- 
stantly decreases, accordingly, its contribution to the 
analogous series circuit resistance constantly 
decreases. The net effect is that as the gap decreases, 
the total reluctance of the magnetic circuit constantly 
decreases. Therefore, the inductance constantly 
increases monotonlcally. For a given change in gap, the 
rate of change of inductance is greatest when the gap is 
the smallest. Accordingly, a system according to a pres- 
ently preferred embodiment has the desirable charac- 
teristic that it is most sensitive to changes in armature 
position when the gap is the smallest, thus enabling the 
most refined control where it is needed the most, i.e., 
when the amriature is close to striking the stator core. 
[0032] The remainder of the inductance computa- 
tion depicted schematically in Figure 5 is designed to 
account for the non-linear way in which flux builds with 
respect to the cun-ent and the gap. The non-linear flux 
characteristic are functions of the air gap and the mag- 
netic permeability of the materials used to fabricate the 
actuator. Because the magnetic permeability of the 
materials will vary depending on the particular alloys 
used, heat treating applied, and other related factors, in 
a preferred embodiment, two independent approxima- 
tions may be applied to account for the air gap and var- 
iable permeability. 

[0033] The first independent approximation is 
termed the "gap factor" approximation. The gap factor 
accounts for the non-linearity of the effect of the gap on 
magnetic flux. This approximation is necessary 
because the flux density is a function of gap size. The 
second independent approximation, accounting for the 
non-linearity of the B-H saturation characteristic, is 



termed the : factor (or "mu" factor) approximation. The 
mu factor approximation accounts for the non-linear 
permeability of ferromagnetic materials. 

5 " Mu" Factor Compensation 

[0034] The magnetic flux density, B, is related to the 
magnetic field intensity, hi. according to the equation 
B = :H = :^:(^H , where Iq is the permeability of space 

10 ( : Q =4 X 1 0 * hen rys/meter ) and measures the effect 
of the magnetic dipole moments of the atoms compris- 
ing the material. The B-H characteristic is a function of 
the magnetic properties of the materials used to fabri- 
cate the actuator. A typical B-H characteristic for ferro- 

15 magnetic materials is depicted in Figure 7. The B-H 
characteristic demonstrates graphically that permeabil- 
ity of ferromagnetic materials varies in a non-linear fash- 
ion as magnetic field strength changes. Referring to 
Figure 7, as magnetomotive force is applied to a mag- 

20 netic circuit, the magnetic flux density increases in a 
non-linear fashion up to the point where the magnetic 
material reaches saturation and the curve begins to 
level off. 

[0035] A table of mu factors for different air gaps 

25 can be constructed as follows. During the time the 
annature 24 is at rest against a pole piece 14, the cur- 
rent may be ramped up and down, taking care to avoid 
allowing the current to drop below the threshold 
required to maintain the amriature 24 in contact with the 

30 pole piece 14. As the current changes, the coil voltage 
may be sampled and, together with the associated cur- 
rent level for each sampled voltage, used to compute a 
table of inductance values associated with each sam- 
pled voltage and current level. From the table of induct- 

35 ance values, a table of mu factors, characteristic of the 
B-H curve of the material used to fabricate the actuator, 
may be readily obtained. The above-described calibra- 
tion process may be performed while the actuator is 
installed and operating in its intended environment. For 

40 example, in the case of an EVT actuator, the calibration 
may be performed while an engine is running while the 
actuator is in a "valve-open" position by varying the cur- 
rent and measuring the corresponding coil voltages, as 
described above. 

45 [0036] The above described mu factor calibration 
may be perfomried on every actuator cycle, or less fre- 
quently, as desired. Once calibrated for a particular 
actuator, the mu factors will typically not change dramat- 
ically from minute-to-minute. However, the mu factors 

50 will tend to vary with temperature and the age of the 
actuator. 

Gap Factor Compensation 

55 [0037] The gap factor accounts for changes in the 
B-H characteristic as the armature moves within the 
actuator As depicted in Figure 7, the shape of the B-H 
curve depends on the air gap of the actuator. As the 
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armature moves within the solenoid, the relative penrie- 
ability of the system changes due to changes in the 
number of lines of magnetic flux coupled through the 
annature. The change in relative permeability in-turn 
changes the B-H characteristic of the system. The gap 
factor approximation accounts for the change In relative 
permeability. The gap factor is not measured directly; 
rather, the gap factor is successively approximated as 
being inversely proportional to the distance between the 
armature and the stator core. 

[0038] The gap factor approximation is founded on 
the principle that when the gap is zero, the full effect of 
the B-H curve is felt by the amnature because permea- 
bility of the solenoid core is maximum. Conversely, 
when the gap is very large, there is only air in the mag- 
netic circuit and the average relative permeability of the 
solenoid core is at a minimum due to the large reluc- 
tance gap with a pemneability of air. As shown in Figure 
7, when the air gap is large, the effect of the B-H curve 
on the armature is minimized. The variation of the B-H 
curve effect between zero-gap (all metal) and a very 
large gap (e.g., an air gap of several centimeters) may is 
approximated in a preferred embodiment as obeying an 
inverse relationship (i.e., a 1/x relationship). 
[0039] The gap factor may be estimated during the 
armature strolce by a succession of approximations as 
follows. A first estimation of inductance L is made, 
assuming ideal gap factors. The estimated value of L 
may then be fed back to estimate the actual (non-ideal) 
gap factor necessary to produce the first estimated 
value of L. The process is repeated to successively 
refine the gap factor until the process converges to zero 
gap under the full effect of the B-H curve. This tech- 
nique offers the benefit of progressively better position 
estimation as the armature 24 approaches the stator 
14. Accordingly, maximum stator control may be 
achieved during the critical period when the arma- 
ture/pole piece gap is on the order of tens of microns 
and the full effect of the B-H curve is realized. 
[0040] Referring to Figure 5, after scaling induct- 
ance, L 70, to units of milli-henrys 72, the inductance 
signal, L, may be compensated 90 by the mu 76 and 
gap 78 factors. After correcting for the gap factor and 
mu factors, inductance, L, is normalized, as depicted in 
88 of Figure 5, to vary preferably between near zero at 
a large gap to a maximum value of 1 .0 at zero gap. The 
maximum inductance may be normalized to any 
number, 1 .0 was chosen in this embodiment for conven- 
ience. The normalization of L accounts for variations in 
absolute inductance that may exist between different 
actuators of like design. Nonnalizing inductance also 
has the benefit of standardizing the range of input sig- 
nals expected by the control system that receives the 
nonnalized inductance as an input For example, the 
actual inductance of a particular actuator may range 
from 10 mH, at maximum gap to 35 mH at zero gap, 
white the actual inductance of a different actuator of like 
design may range from 12 mH at maximum gap to 40 



mH at zero gap. Normalizing the inductance allows for 
automatic calibration between actuators of different 
absolute inductance and simplifies the control loop 
design for a standard range of inputs. 

5 

Velocity State Variable Estimation 

[0041] As with the armature position estimation, 
which was derived by dynamically measuring induct- 

10 ance, as described above, the velocity state variable 
may be estimated by measuring rate of change of 
inductance. The "brute force" approach of differentiating 
the position signal to obtain the armature velocity does 
not generally achieve satisfactory results because 

15 minor "noisy" perturbations inherent in the position sig- 
nal will have very large derivatives, and hence, will pro- 
duce a corrupt velocity signal. Accordingly, armature 
velocity must be measured by an atternative method. 
[0042] In accordance with a presently preferred 

20 embodiment, the armature velocity may be approxi- 
mated by investigating the integral-derivative relation- 
ship between rate of change of magnetic flux, d<t>/dt, 
and magnetic flux, 4>. and recognizing that dL/dt is pro- 
portional to armature velocity as follows. As described 

25 above, position may be estimated by mapping induct- 
ance. L. to position, where L is, in-turn, determined by 
dividing flux, <t> by coil current. I. in accordance with the 
expression <t> = LI. in like fashion, armature velocity may 
be directly estimated from the d4>/dt signal, as calcu- 

30 lated at 62 in Figure 5. Because d<t>/dt is a relatively 
uncorrupted "clean signal" it may be used as a suffi- 
ciently precise estimate of armature velocity to enable a 
soft landing of the amnature against the stator core. 
[0043] The velocity approximation derivation is as 

35 follows: Given the basic relationship of 4> = LI, d(t>/dt 
may be approximated as I dL/dt, where "1" is a real time 
measured value of the instantaneous current magni- 
tude, therefore dl/dt does not have to be considered. 
Accordingly dL/dt is approximately equal to (d<I>/dt)/l. In 

40 a preferred embodiment, dL/dt may be scaled by the 
same mu factor and gap factor used to scale L. The 
result of scaling dL/dt by the mu and gap factors is 
labeled "du/dt" in the present disclosure ("du/dt" is a 
"dummy variable" representing a rate term) and may be 

45 used to approximate armature velocity. 

[0044] Referring again to Figure 5, the above 
method may be implemented by dividing dcl>/dt, the out- 
put of 62 (d<t>/dt) by coil cun-ent, I, at 74. The resulting 
approximated value of dL/dt may then be compensated 

50 by the mu 76 and gap 78 factors, and scaled by a con- 
stant 82 to produce a rate temn. du/dt 84. corresponding 
to armature velocity. Accordingly, the outputs of the sys- 
tem depicted in Figure 5 are normalized inductance, L 
86, (the position estimation term) and rate of change of 

55 inductance, du/dt, (the velocity estimation term). 

[0045] Accordingly, in a presently preferred embod- 
iment, the inductance L may be detenmined by measur- 
ing the magnetic flux. The rate of change of inductance 
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may be estimated as being proportional to the rate of 
change of flux. The resulting state variables constitute 
the inputs to a control system for modulating coil cur- 
rent, and hence controlling armature velocity. A signifi- 
cant benefit of the above-described system for 
dynamically estimating the actuator state variables of 
position, velocity, and magnetic flux density is that there 
is no differentiation required to obtain rate information. 
For the reasons given above, it is highly desirable to 
avoid differentiation of non-ideal signals. 
[0046] Another feature of the above-described 
approach is the ability to compensate inductance L for 
variations in the B-H characteristic due to changing per- 
meability as the armature moves within the solenoid. 
Normalizing L and obtaining armature rate information 
from rate of change of flux, dct>/dt, also contributes to the 
simplicity of the actual implementation. 
[0047] The function of the mu factor can best be 
appreciated with reference to Figure 8. Figure 8 depicts 
typical data obtained during the autocalibration of the B- 
H cun/e and mu factor table loading. As the current is 
ramped up and down, the inductance changes in 
inverse proportion to the current through the coil. 
Accordingly, as the current decreases, the inductance 
increases. 

[0048] Waveform 110 in Figure 10 is the integral 
rate signal (the dL/dt signal in a preferred embodiment) 
and waveform 1 1 2 is the estimated inductance L. Note 
that the shapes of the curves are very similar, thus vali- 
dating empirically the simplifying assumptions that 
d4>/dt may be approximated as I dLVdt, and dl/dt is neg- 
ligible. These assumptions greatly reduce the complex- 
ity of the implementation hardware and/or software. 

Estimation Of Closed Loop Controller Set Points (Run- 
ning Set Points) 

[0049] To this point we have set forth sense rless 
methods for obtaining the state variables of magnetic 
flux, armature position and armature velocity. It remains 
to be described how the state variables are used by a 
control system to control armature velocity and gener- 
ate a soft landing against a stator core. The control sys- 
tem must receive as inputs armature position and 
velocity information to achieve a soft landing. In addi- 
tion, as the armature approaches the stator core, 
greater precision and accuracy are required in the posi- 
tion and velocity estimates. Closed loop controller set 
points provide continuously updated target positions 
and velocities during the armature stroke. 
[0050] Several basic principles governing control 
system design have become apparent during experi- 
mental testing. First, the control system should not start 
attempting to control armature velocity until the arma- 
ture moves close enough to the stator core such that 
there is sufficient flux passing through the armature to 
exert significant control over the armature by changing 
the coil current. Stated another way, there must be suf- 



ficient magnetic energy in the working gap before the 
control system can exert control over the armature. As a 
rule of thumb, the armature should be close enough to 
the stator core that the amount of magnetic flux closed 
5 through the core is at least equal to the amount of flux 
that escapes the core. Attempting to exert control over 
the armature before sufficient magnetic flux has been 
closed through the core will result in ineffective control, 
large coil current and associated power dissipation in 
the form of heat. 

[0051] The reluctance path of the actuator corre- 
sponds to the armature air gap. As explained above, 
reluctance is analogous to resistance in dc-resistive-cir- 
cuit analysis and is defined as the ratio of the magneto- 

15 motive force to the total flux. When the air gap is large, 
the reluctance is great and a large portion of magnetic 
flux will leak away and not pass across the air gap 
where it is needed to control the force on the armature. 
Accordingly, it is ineffective to close the control loop on 

20 the system until the air gap is sufficiently small (i.e., the 
armature is close to the stator core) to keep flux from 
leaking away from the air gap. 

[0052] When the armature is sufficiently close to 
the stator core for the system to exert control over the 

25 armature by varying the magnetic flux in the circuit, the 
control circuitry "closes the loop" and begins controlling 
the armature velocity. Once the armature is placed 
under closed loop control, running set points are deter- 
mined con^esponding to intermediate armature position 

30 and velocity targets during the armature stroke. The 
term "running set point" refers to a control system target 
for position or velocity that changes dynamically during 
the armature stroke. As the armature moves towards 
the stator core under closed-loop control, the set points 

35 for position and velocity are dynamically updated until 
the armature lands on the stator core (i.e., zero veloc- 
ity). Running set points can be thought of as defining a 
near-optimal armature position and velocity trajectory 
sufficient to achieve a soft landing of the armature 

40 against the stator core. 

[0053] Figure 1 1 depicts the normalized inductance 
and rate of change of inductance that may be empiri- 
cally determined as the optimal values for the running 
set points. The loop closed at 1 14 with initial set points 

45 116 and 117. Under multivariate closed loop control, 
armature velocity 118 decreases as the set points are 
updated 120 and 121. As the armature continues to 
move, the system follows the updated set points 120 
and 121 until the armature lands at near-zero velocity 

50 122. 

[0054] Figure 1 1 also demonstrates that in an alter- 
native preferred embodiment, the control loop logic may 
use inductance and rate of change of inductance 
directly as the state variables for controlling the system. 
55 In this embodiment, a reduction in hardware complexity 
is achieved because there is no need to mathematically 
convert inductance and rate of change of inductance 
into respective position and velocity terms for input to 
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the control loop logic. Set point traces 1 20 and 1 21 in 
Figure 11 demonstrate application of this method; 
rather than position and velocity, the traces represent 
the inductance and rate of change of inductance inputs 
to the control loop logic. After It was experimentally 
demonstrated that armature position and velocity could 
be accurately estimated from inductance and rate of 
change of inductance, it was further determined that 
amriature velocity could be placed under multivariate 
control based directly on inductance and rate of change 
of inductance. Accordingly, the set points 120 and 121 
of Figure 1 1 are actually in units of inductance and rate 
of change of inductance, rather than position and veloc- 
ity. 

[0055] Figure 12 is a block diagram demonstrating 
how running set points may be determined using nor- 
malized inductance and rate of change of inductance as 
inputs. The actual set point target values for inductance 
and rate of change of inductance are determined empir- 
ically and adjusted over the entire armature stroke to 
achieve an ideal soft landing of the armature against the 
statorcore. The ideal set point values are stored in look- 
up tables, represented as 130 and 132. Note the "posi- 
tion" and "velocity" set point tables in Figure 12, 130 and 
132, respectively, may also correspond with inductance 
and rate of change of inductance in accordance with an 
alternative preferred embodiment described above. The 
set points may be empirically detemriined by adjusting 
an actuator for a perfect soft landing and recording the 
ideal trajectory of normalized inductance and rate of 
change of inductance. 

[0056] The set points represent the Ideal position 
and velocity (or, in a preferred embodiment, inductance 
and rate of change of inductance) of the armature at 
every point in the stroke. As depicted in Figure 12, dur- 
ing operation, the actual normalized inductance 131 (or 
position in an alternative embodiment) is subtracted 1 34 
from the appropriate set point corresponding to induct- 
ance (or position in an altemative embodiment), yielding 
a proportional error 136. In similar fashion, the rate of 
change of inductance 133 (or velocity in an alternative 
embodiment) is subtracted 138 from the appropriate set 
point corresponding to rate of change of inductance (or* 
velocity in an alternative embodiment), yielding a corre- 
sponding rate error 140. The proportional error 136 and 
rate error 140 at multiple instants of time may then be 
applied as inputs to the control system logic. 

Control System Logic 

[0057] Figure 13 is a comparison of measured 
inductance 1 42 with ideal inductance 1 44 in accordance 
with a presently preferred embodiment. In this example, 
a conventional PID (proportional, integral, derivative) 
servo was used in this example to demonstrate the fea- 
sibility of tracking the ideal set point values of induct- 
ance. The control loop used the proportional error 
signal 143 as a feedback input. It may also be observed 



that under closed loop control, the PID controller varied 
the current based on the error signal to force the meas- 
ured inductance signal 142 track with the ideal set 
points for inductance 144. 

5 [0058] Figure 15 demonstrates that a soft landing 
was achieved in accordance with the above-described 
methods, using the PID controller system described 
above in reference to Figure 14. In this example, a 
dSPACE, Inc. 1102 commercial DSP microprocessor 

10 controller board, with a Texas Instruments TMS320 
DSP was used. However, any conventional DSP or ana- 
log controller may be substituted. As shown in Figure 
15. the velocity of the armature in region 146, as the 
armature approaches the state r core, is sharply 

15 reduced, thus enabling a soft landing. 

[0059] Any of several known multivariate control 
algorithms may be applied for closing the control loop 
based on the proportional error 136 and rate en'or 140. 
In a preferred embodiment, the control system is a fuzzy 

20 logic controller. In an alternative preferred embodiment, 
the control system is a state feedback system. 
[0060] While the present invention has been dis- 
closed with reference to certain preferred embodiments, 
numerous modifications, alterations, and changes to the 

25 described embodiments are possible without departing 
from the sphere and scope of the present invention, as 
defined in the appended claims. Accordingly, it is 
intended that the present invention not be limited to the 
described embodiments, but have the full scope defined 

30 by the language of the following claims, and equivalents 
thereof. 

Claims 

35 1. A sensorless method of controlling the landing 
velocity of an armature in an electromagnetic actu- 
ator, comprising the steps of: 

providing an electromagnetic actuator having a 

40 coil; 

measuring the inductance of the coil in real- 
time as the armature moves within the coil; 
compensating the measured inductance for 
non-linear permeability and magnetization 

45 effects; and 

providing the measured inductance to a control 
system for modulating a current delivered to 
the actuator. 

50 2. The method of claim 1 , further comprising the step 
of normalizing the measured inductance at zero 

gap- 

3. The method of claim 2, further comprising the steps 
55 of: 

estimating the rate of change of inductance of 
the coil in real-time as the armature moves 
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within the actuator; 
compensating the estimated rate of change of 
inductance for non-linear permeability and 
magnetization effects; and 
providing the compensated rate of change of 
inductance to a control systenn for modulating a 
current delivered to the actuator. 

4. The method of claim 3, wherein the rate of change 
of inductance is determined without differentiating 
the inductance signal. 

5. The method of claim 4, further comprising the step 
of capturing the B-H magnetization characteristics 
of the actuator during an armature stroke. 

6. The method of claim 5, wherein the step of captur- 
ing the B-H magnetization characteristics of the 
actuator during an armature stroke further includes: 

maintaining the armature in contact with a pole 
piece; 

driving a time-varying current through the coil; 
sampling the voltages associated with a plural- 
ity of current levels; 

computing inductance values associated with 
each sampled voltage and current level; and 
computing mu factors for each inductance 
value. 

7. The method of claim 5, wherein the inductance cor- 
responds to an armature position estimation and 
the rate of change of inductance corresponds to an 
armature velocity estimation. 

8. The method of claim 5, further comprising the step 
of measuring the coil resistance of the actuator 
while the armature is in a rest position against a 
stator core. 

9. The method of claim 8, wherein the step of measur- 
ing the coil resistance of the actuator while the 
armature is in a rest position against a stator core 
further includes: 

driving the coil with a steady-state current; 
measuring the coil voltage necessary to main- 
tain the steady-state current through the coil; 
and 

dividing the measured voltage by the steady- 
state current to calculate coil resistance. 

10. The method of claim 9, wherein the control system 
is a fuzzy logic control system. 

11. The method of claim 9, wherein the control system 
is a full state feedback control system. 



18 

12. The method of claim 9. wherein the control system 
is a PID control system. ' 

13. The method of claim 9, wherein the electromechan- 
5 leal actuator is operatively attached to a fuel injec- 
tor. 

14. The method of claim 13, wherein the fuel injector is 
a direct injection fuel injector. 

10 

1 5. The method of claim 9, wherein the electromechan- 
ical actuator is an EVT actuator. 

16. The method of claim 9, wherein the control system 
15 comprises a microprocessor. 

17. The method of claim 9, wherein the control system 
comprises a digital logic circuit. 

20 18. The method of claim 9, wherein the control system 
comprises an analog circuit. 

19. A method of controlling the velocity of an armature 
in an electromagnetic actuator as the amiature 

25 moves from a first position towards a second posi- 
tion, the electromagnetic actuator including a coil 
and a core at the second position, the coil conduct- 
ing a current and generating a magnetic force to 
cause the annature to move towards and land at 

30 the second position, and a spring structure acting 
on the armature to bias the armature from the sec- 
ond position, the method comprising the steps of: 

measuring the inductance of the coil as the 
35 armature moves within the actuator; 

compensating the measured inductance for 
non-linear permeability and magnetization 
effects; and 

providing the measured inductance to a control 
40 system for modulating a current delivered to 

the actuator. 

20. The method of controlling velocity of an armature in 
an electromagnetic actuator according to claim 1 9, 

45 further comprising the step of normalizing the 
measured inductance at zero gap. 

21. The method of controlling velocity of an armature in 
an electromagnetic actuator according to claim 20. 

50 further comprising the steps of: 

estimating the rate of change of inductance of 
the coil as the armature moves within the actu- 
ator; 

55 compensating the estimated rate of change of 

inductance for non-linear permeability and 
magnetization effects; and 
providing the compensated rate of change of 
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inductance to a control system for modulating a . 29. The method of claim 27, wherein the control system 
cun-ent delivered to the actuator. is a full state feedback control system. 



22. The method of controlling velocity of an amnature in 
an electromagnetic actuator according to claim 21, 
wherein the rate of change of inductance is deter- 
mined without differentiating the inductance signal. 

23. The method of controlling velocity of an amiature in 
an electromagnetic actuator according to claim 22, 
further comprising the step of capturing the B-H 
magnetization characteristics of the actuator during 
an armature stroke. 

24. The method of controlling velocity of an armature in 
an electromagnetic actuator according to claim 23, 
wherein the step of capturing the B-H magnetiza- 
tion characteristics of the actuator during an arma- 
ture stroke further includes: 

maintaining the armature in contact with a pole 
piece; 

driving a time-varying current through the coil; 
sampling the voltages associated with a plural- 
ity of current levels; 

computing inductance values associated with 
each sampled voltage and cun^ent level; and 
computing mu factors for each inductance 
value. 

25. The method of controlling velocity of an annature in 
an electromagnetic actuator according to claim 23, 
wherein the inductance corresponds to an anna- 
ture position estimation and the rate of change of 
inductance corresponds to an annature velocity 
estimation. 

26. The method of controlling velocity of an annature in 
an electromagnetic actuator according to claim 23, 
further comprising the step of measuring the coil 
resistance of the actuator while the armature is in a 
rest position against a stator core. 

27. The method of controlling velocity of an annature in 
an electromagnetic actuator according to claim 26, 
wherein the step of measuring the coil resistance of 
the actuator while the armature is in a rest position 
against a stator core further includes: 

driving the coil with a steady-state current; 
measuring the coil voltage necessary to main- 
tain the steady-state current through the coil; 
and 

dividing the measured voltage by the steady- 
state current to calculate coil resistance. 

28. The method of claim 27, wherein the control system 
is a fuzzy logic control system. 



30. The method of claim 27, wherein the control system 
5 is a PID control system. 

31. The method of claim 27, wherein the electrome- 
chanical actuator is operatively attached to a fuel 
injector. 

10 

32. The method of claim 31 , wherein the fuel Injector Is 
a direct injection fuel injector. 

33. The method of claim 27. wherein the electrome- 
15 chanical actuator is an EVT actuator. 

34. The method of claim 27, wherein the control system 
comprises a microprocessor. 

20 35. The method of claim 27, wherein the control system 
comprises a digital logic circuit. 

36. The method of claim 27, wherein the control system 
comprises an analog circuit. 

25 

37. An apparatus for controlling velocity of an armature 
in an electromagnetic actuator as the armature 
moves from a first position towards a second posi- 
tion, the electromagnetic actuator including a coil 

30 and a core at the second position, the coil conduct- 
ing a current and generating a magnetic force to 
cause the annature to move towards and land at 
the second position, and a spring structure acting 
on the armature to bias the armature from the sec- 

35 ond position, the apparatus comprising: 

a means for estimating the rate of change of 
inductance of the coil as the armature moves 
within the actuator; 
40 a means for compensating the estimated rate 

of change of inductance for non-linear permea- 
bility and magnetization effects; 
a means for normalizing the measured induct- 
ance at zero gap. 
45 a means for estimating the rate of change of 

inductance of the coil in real-time as the arma- 
ture moves within the actuator; 
a means for compensating the estimated rate 
of change of inductance for non-linear permea- 
se bility and magnetization effects. 
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electromagnetic solenoid 



(57) An improved nnethod for controlling the landing 
velocity of an armature in an electromechanical actua- 
tor, such as a fuel injector, fuel pressure regulator, or 
engine valve actuator is provided. The position and ve- 
locity of an annature during a stroke is dynamically es- 
timated by calculating the inductance and rate of change 
of inductance of the actuator coil in real-time as the ar- 
mature moves through its stroke, compensating for non- 
linear permeability and magnetization effects due to 
changing gap, temperature, magnetic material proper- 
ties or magnetic architecture, normalizing the calculated 
inductance value at the end of a stroke (zero gap), and 
mapping the value of normalized Inductance to corre- 
spond to an armature position by an algebraic transfor- 
mation. Inductance may be used directly as a position 
variable without mapping it to units of position. Rate of 
change of inductance may be used as a rate variable 
without mapping it to units of velocity. 
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